Experimental investigation of the initial regime in fingering electrodeposition: 
dispersion relation and velocity measurements 
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Recently a fingering morphology, resembling the hydrodynamic Safi'man- Taylor instability, was 
identified in the quasi-two-dimensional electrodeposition of copper. We present here measurements 
of the dispersion relation of the growing front. The instability is accompanied by gravity-driven 
convection rolls at the electrodes, which are examined using particle image velocimetry. While at 
the anode the theory presented by Chazalviel et al. describes the convection roll, the fiow field 
at the cathode is more complicated because of the growing deposit. In particular, the analysis of 
the orientation of the velocity vectors reveals some lag of the development of the convection roll 
compared to the finger envelope. 

PACS numbers: 81.15Pq, 47.54, 89.75Kd 



I. INTRODUCTION 

It is sometimes believed, that all interesting phenom- 
ena in the universe happen at interfaces Q. Following 
this line of thought, we believe that the study of the 
dynamics of interfaces provides a key for understand- 
ing generic features of non-equilibrium phenomena. The 
electrochemical deposition of metals from aqueous solu- 
tions in quasi-two-dimensional geometries is an easily ac- 
cessible growth phenomenon of such an interface. The 
emerging structures show a broad variety of growth pat- 
terns including fractals, seaweed or dendrites. For a re- 
cent review see and references therein. 

The focus of this paper is on the electrodeposition of 
finger deposits |3|: after the addition of a small amount 
of an inert electrolyte like sodium sulphate to a cop- 
per sulphate solution, the morphology of copper deposits 
changes from a typical fractal or dense-branched red cop- 
per structure to some fine-meshed texture with a finger- 
like envelope. Figure gives an example of the early stage 
of a deposit formed under these circumstances. The un- 
derlying mechanism is believed to be qualitatively under- 
stood. The increase of the electric conductivity enables 
alternative reaction paths like the reduction of H2O. The 
resulting increase of the pH value triggers the formation 



of a copper hydroxide gel (Cum(OH)„ 



(2m-n) + N 



in front 



of the advancing deposit |J, |5[. When considering, that 
the fluid between the copper filaments contains no gel, 
the ensuing situation resembles the Saffman- Taylor in- 
stability, where a more viscous fluid is pushed by a less 



viscous one and their interface develops the same type of 
fingering (see Ref. D for a recent survey). 

The Saffman- Taylor instability is strongly influenced 
by the surface tension of the interface. In this paper we 
use that idea to measure the strength of an effective sur- 
face tension associated to the hydrogel- water interface by 
analyzing the dispersion relation. It should be remarked, 
that the nature of surface tension between miscible fluids 
is still an active area of research 0| . 

Another necessary ingredient for the occurrence of fin- 
gers are density-driven convective currents in front of the 
growing deposit: If convection is suppressed by turning 
the electrodeposition cell in a vertical configuration, fin- 
gers are not longer formed ||. For that reason it ap- 
pears to be essential to understand the nature of the con- 
vection field in our experiment, which we examine using 
particle image velocimetry (PIV). 

The organization of the paper is the following: In Sec- 
tion |^ we introduce the experimental setups. Section III 
is devoted to the dispersion relations, with [II A cover- 



ing some technical aspects and [II B presenting the mea- 
sured dispersion relations and analyzing the results for 
a textured electrode. In Section [V we discuss the PIV 
measurements: [VA is devoted to the anode while [VB 
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summarizes our results for the cathode. Finally, Section 
^ contains our conclusions. 



II. EXPERIMENTAL SETUPS 

The electrodeposition is performed in a cell with two 
glass plates of 8 x 8 cm^ area as side- walls. Two par- 
allel copper wires (99.9 %, Goodfellow) separated by a 
distance of 4 cm serve as electrodes and spacers. Their 
diameter d ranges between 125 /xm and 300 /im. Figure 




FIG. 1: The Finger morphology 403 s after the start of the 
experiment. The solid black line at the bottom of the image 
is the cathode. The width of the image corresponds to 22.1 
mm. The deposit was grown in a cell of 250 /im thickness, 
the applied potential was 15 V. The electrolyte contained 50 
mM CuS04 and 4 mM Na2S04. 
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TABLE L Image acquisition systems used in the experimental 
setups. 
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^ a) shows a sketch of this setup. We use a coordinate 
system where the x axis is parallel to the electrode, the y 
axis points from the cathode to the anode and the z axis 
is perpendicular to the glass plates. 

The space between the electrodes is filled with an aque- 
ous solution of 50 mM CUSO4 and 4 mM Na2S04 for 
the measurements of the dispersion relations and 50 mM 
CUSO4 and 7 mM Na2S04 for the PIV experiments. All 
solutions are prepared from Merck p. a. chemicals in non 
deaerated ultrapure H2O. 

All measurements are performed with constant (within 
0.4 %) potential between the electrodes ranging between 
12 V and 19 V. The average current density is below 35 
mA / cm^ . 

The two targets of our investigation require two differ- 
ent ways of illumination (sketched in Fig. p|(b) and (c)) 
and image acquisition (summarized in tableffl). Since the 
measurements of the dispersion relation demand a high 
spatial resolution we used a Kodak Megaplus 6.3i CCD- 
camera with 3070 x 2048 pixel mounted on a Nikkor SLR 
macro lens with 105 mm focal length and a spacer ring. 
To take full advantage of the spatial resolution of 7.9 /im 
per pixel it was necessary to employ a Kohler illumina- 
tion Q using filtered light with a wave length of 405 nm 
from a tungsten lamp. Images are taken in intervals Ataq 
of 5 s and are directly transferred with a frame-grabber 
card to the hard disk of a PC. 

To visualize the velocity field inside the cell, we added 



FIG. 2: (a) Electrodeposition cell with finger deposit and our 
choice of the coordinate system, (b) Side view of the illumi- 
nation with shining-through light used for the measurements 
of the dispersion relation. Image |l] was taken that way. (c) 
Side view of the dark field microscopy used for the PIV. Im- 
age ^ gives an example of this illumination technique. The 
ovals inside the cell indicate the convection rolls. 



latex tracer particles to the electrolyte. We used parti- 
cles with 0.3 fim diameter, which stay suspended due to 
Brownian motion. Since we cannot resolve these particles 
with our optical system, we used dark-field microscopy: 
only light scattered from objects inside the cell falls into 
the lens. Imaged gives an example, the white area at the 
bottom represents the growing deposit, the points above 
correspond to tracer particles. 

We did not observe electro-osmosis as reported in [||, 
but the particles show some tendency to coagulate and 
settle to the bottom plate. This problem is handled in 
later stages of the image processing. 

Images were acquired using an Olympus SZH stereo 
microscope and a Sony XC 77RR CE CCD-camera with 
512 X 512 pixels which resulted in a spatial resolution 
of 17 /zm per pixel. Ataq was 2 s and images were also 
directly transferred to a PC. 




FIG. 3: Part of an image taken in dark-field microscopy. Only 
light scattered from the deposit or the tracer particles is vis- 
ible. The width corresponds to 3.4 mm. (Cell thickness: 300 
fim, applied potential: 12 V, solution: 50 mM CuS04 and 7 
mM Na2S04.) 



III. DISPERSION RELATION 

To characterize the instability of a pattern forming sys- 
tems there is a quite common method: starting with the 
uniform system, one adds a small sinusoidal perturba- 
tion of wavenumber k and amplitude Aq and investigates 
its temporal evolution. As long as the system is in the 
linear regime the perturbation will grow or shrink expo- 
nentially: 

A(fc,t) = Ao e''('=)* (1) 

The dependency of the growth rate a on the wavenumber 
k is called dispersion relation. 

In the context of electrodeposition dispersion relations 
have been measured for the initial phase of compact pC| ] 
and ramified |0 growth and calculated to explain the 
stability of the dense radial morphology |jl^ . 

A. Image processing 

In order to measure the dispersion relation the first 
step is to track down the temporal evolution of the front: 
in the image taken at time t we identify for each column 
X the height h{x,t) of the deposit. This is done in two 
steps: First we search for the pair of pixels with the 
highest gray value gradient in each column. Second we 
perform a sub-pixel interpolation by calculating the point 
of inflexion between these two pixels. By evaluating a 
stationary edge, we could demonstrate that the error of 



3 

this algorithm is smaller than 0.1 pixel. However the 
random fine structure of the deposit acts like the addition 
of shot noise to h{x,t). We try to mitigate this effect by 
applying a median filter with 7 pixel width ||l^ . Figure 
^ shows h{x, t) for 4 various times of an experiment. 
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FIG. 4: Result h{x, t) of the front tracking. The curve at 
the bottom corresponds to the initial position of the cathode. 
Above the front at 156 s, 259 s, 362 s and 466 s after the 
beginning of the experiment. The experiment is the same as 
the one presented in Fig. |l[ The part of length w between the 
dashed lines was chosen by the cut-out algorithm for further 
processing. 



Our next aim is to do a Fourier decomposition of ft, (a; , t) 
to examine the dynamics of the amplitudes A{t^m) of 
the individual modes to. Because our initial conditions 
are random noise, which can be understood as a super- 
position of various modes with different wavelengths A, 
we encounter a leakage problem: We observe our sys- 
tem in a window of width wq. The ratio wq/X of the 
dominant modes is almost always fractional. This has 
a very unpleasant consequence for the Fourier transfor- 
mation, which computes the amplitudes for integer mode 
numbers: power of the dominant modes is transfered to 
other less excited modes and spoils the measurement of 
the growth rates there. 

A common answer to this issue is the usage of a win- 
dowing function |p^ , which does not eliminate leakage, 
but restricts it to neighboring mode numbers. The price 
to be paid is, that this leakage will now occur even if wq 
is a common multiple of all wavelengths contained in the 
initial signal. 

In this paper we use a different approach, which is 
illustrated in Fig. |[ we cut out a part of h{x, t) of width 
w and offset o and use only this part for the subsequent 
analysis, o and w, which are constant for the whole run 
of the experiment, are chosen such that the left and the 
right end of the cut-out have the same height: 

h{o,t) - h{o + w,t) ^0 (2) 



4 



and the same slope: 

dh{o,t) 
dx 



dh{o + w, t) 
dx 



1000 



= 



(3) 



Fulfilling Eq. || and Eq. ^ for all times is for all practical 
purposes identical to the statement that w is a common 
multiple of all wavelengths contained in h{x,t), that is 
the condition under which no leakage occurs. 

In practice, Eq. |^ and Eq. ^ can be satisfied only ap- 
proximately. The algorithm followed to select o and w 
minimizes the sum over the height differences according 
to Eq. H while it assures, that the sum over the slope 
differences according to Eq. || does not exceed a thresh- 
old. The average height difference obtained in that way 
is < 3 pixels while the slope difference threshold was 1 
pixel/pixel, both values correspond to the noise level of 
h(x,t). The final step is a standard Fourier decomposi- 
tion of the cut-out parts of h(x,t). 



B. Results 

For each Fourier amplitude we tried to fit the temporal 
evolution with the exponential growth law given by Eq. |^. 
Fig. 1^ gives an example for two modes, corresponding to 
the experiment displayed in Fig. ^. The time interval 
for the exponential fits is indicated by the solid symbols. 
The start time is set by the time the hydrogel layer needs 
to build up; before that point no instability or destabi- 
lization of the growing interface is clearly evidenced. In 
the experiment shown in Fig. ^ this point corresponds to 
an elapsed time of about 150 s. However some modes 
need longer until they are grown to an amplitude which 
can be distinguished from the measurement background 
noise. At the other end the fit is limited by the onset of 
deviations from the exponential growth law due to non- 
linear effects becoming important. 

To test the dependence of the dispersion relation on 
the applied potential, we performed measurements with 
15 V and 19 V using electrodes of 250 /im diameter. The 
results are displayed in Fig ^, each is averaged over three 
experiments. Both dispersion relations show a limited 
band of positive amplitude growth rates with k between 
zero and fccrit and a wave number /cmax where the growth 
rate is maximal. 

The increase of U from 15 V to 19 V is accompanied 
by an increase of the average growth velocity v from 13.3 
/xm/s to 15.9 /xm/s. The distinct shift of fcmax to higher 
wave numbers is in conformity with measurements of the 
number of incipient fingers as a function of v reported in 
Q . The decrease of /ccrit could originate from a change of 
the physicochemical properties of the copper hydrogel. 

Negative "growth rates" can only be measured if the 
initial amplitude of the corresponding mode is strong 
enough. To enforce this we prepared a textured electrode 
which stimulates the initial growth at a wavenumber of 
62.8 cm~^. It consists of a synthetic substrate of 120 ^m 
height with a 35 /im copper plating. The copper layer 



100 



10 





1 1 1 








- 








□ 


■F ' 






/m 




q W/ 


a* ■ 




100 200 


300 400 500 
Elapsed Time (s) 


600 



FIG. 5: Exponential growth of the Fourier modes. The wave 
numbers are 6.9 cm~^ (circles) and 13.8 cm~^ (squares). Only 
the data points represented by filled symbols were included in 
the fits with Eq. |l]. The experiment is the same as in Fig. ^ 
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FIG. 6: Dependence of the dispersion relation on the applied 
potential. Filled circles correspond to experiments performed 
applying 15 V, open circles applying 19 V. Cell thickness was 
250 /im in both cases. All data sets are averaged over three 
experiments, error bars give the standard deviation of the 
mean value. 



was etched to derive a comb-like structure with copper 
stripes of width 0.75 mm and spacings of 0.25 mm. Is- 
lands of growing deposit evolve at the tips of the copper 
stripes and amalgamate after some time. Fig. ^ gives an 
example. The temporal evolution of two Fourier modes 
of this experiment is shown in Fig. |[ 

The dispersion relation of the textured electrode is dis- 
played in Fig ^ together with results for cell thickness of 
125 /im and 250 /im. Within the scope of our experimen- 
tal errors no influence of d is observable. In principle a 
decrease of d is accompanied by a decrease of the surface 
tension forces, which should result in a higher fccrit- How- 
ever as discussed in chapter [V convection will also be 



FIG. 7; Growing deposit at the textured electrode. Images 
were taken (a) 51 s, (b) 82 s, (c) 142 s, (d) 202s after the 
beginning of the experiment. The width of the images corre- 
sponds to 15.8 mm. (Applied potential: 15 V, cell thickness 
155 ^m, solution: 50 mM CuS04 and 7 mM Na2S04.) 
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FIG. 8: Temporal evolution of the Fourier amplitudes mea- 
sured with the textured electrode. The data belong to the ex- 
periment displayed in Fig. W. The wavenumbers are 3.5 cm~^ 
(squares) and 62.6 cm~^ (circles). Data points represented 
by filled symbols were included in the fits with Eq. hi 



significantly reduced. This will result in a steeper tran- 
sition between the hydrogel and the electrolyte, which 
increases the surface tension 7 ||l5|. Presumably these 
two effects cancel each other. 

While we are not claiming that the Saffman- Taylor dis- 
persion relation: 

Vgel + m-iO 12(?7gei + ??H20) 

gives a full explanation of this fingering phenomenon, we 
do believe that it captures the essentials of the physical 
mechanism, especially the damping effects due to an ef- 
fective surface tension. Therefore we performed a fit of 
Eq. H to the dispersion relation of the textured electrode, 
which is shown in Fig. H. As results we find a viscosity of 
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FIG. 9: Dependence of the dispersion relation on the cell 
thickness. Open circles were measured in cells with d — 250 
fim, open squares in cells with d = 125 ^m. The closed circles 
give the results of the textured electrode (d = 155 fim). The 
sohd line is a fit of Eq. ^ to the growth rates of the textured 
electrode. Applied potential was 15 V in all cases. All data 
sets are averaged over three experiments, error bars give the 
standard deviation of the mean value. 



the hydrogel rjgei of 2.4 lO^'^ kg/ms which is about twice 
the viscosity of water ?7h2 0- The effective surface tension 
7 turns out to be 3.5 10"'' N/m. This is about five mag- 
nitudes lower than the surface tension at the water-air 
interface, which is presumably due to the fact, that the 
gel and the water are miscible fluids. 



IV. VELOCITY MEASUREMENTS 

Elcctrodeposition is often accompanied b y b uoyancy- 
driven convection rolls [|[00^0|o[|ll||. 
The driving force for the convection are the concentra- 
tion changes at the electrodes: at the anode the ion con- 
centration and therefore the density of the electrolyte 
increases. While it descends, lighter bulk solution flows 
in and a convection roll as sketched in Fig || (c) starts to 
grow. 

To visualize the growing deposit in the x-y plane, we 
observe the cells from above. As apparent from Fig. || (c) 
the plane defined by the convection roll is the y-z plane. 
This results in the uncomfortable situation, that the ob- 
served tracer particles move simultaneously towards and 
away from the electrode, which obviates the use of stan- 
dard correlation techniques for the PIV. 

We therefore developed a software package capable to 
keep track of the motion of individual particles. This is 
done in two steps: First we identify the particles in each 
image and insert their center of mass into a database. 
Then we construct contiguous histories for individual 
particles using five consecutive time-steps. In this way 
we are able to measure the Wx and Vy components of the 
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flow. The software is published under the GNU Public 
License and can be downloaded ft-om a URL IG4|. 



Anode Results 



with 



ki ^ 13.1 



dp 



Zc F (^a + Mc) V 



1/3 



1^ 



(6) 



All velocity measurements were performed in cells with 
thickness d = 300 fj,m applying a potential of 12 V. At the 
anode a convection roll is clearly visible: while having no 
relevant velocity component in the x-direction, the tracer 
particles in a distinct zone move towards or away from the 
electrode. Fig. gives the Vy components of all particles 
detected in one image as a function of their distance y to 
the anode. 



and 
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FIG. 10: Particle velocities at the anode after 888 s. Particles 
with positive Vy are moving away from, with negative Vy to- 
wards the anode. Only particles lying outside the solid lines 
were considered for the calculation of an average velocity Wavg . 
The dashed line is a fit of Eq. |l^ to Wavg(3/,t). The applied 
potential is 12 V, the cell thickness 300 /im. 



Due to the big depth of focus of our optical system, 
we observe particles in all heights z of the cell simultane- 
ously. As the particle velocity is a function of z, we find 
for a given y all velocities between ±Vmax{y,t). 



1. Theory 

It is known that some time after the start of an exper- 
iment vertical diffusion starts to smear out the concen- 
tration differences between the flows to and away from 
the electrode. Chazalviel and coworkers proposed 
a two-dimensional description for this diffusion-hindered 
spreading (DHS) regime. The velocity component per- 
pendicular to the electrode Vy should obey [p5[: 



Vy{y,Z,t)^ ki 



1 - 



y 



k2 Vt 



_^f y^ 

12 V fc2 Vt 



- z- 



(5) 



fc, = 0.222 



1/3 



(7) 



i denotes the current density (25 ± 3 mA/cm^), ^ 
is the dependency of the density on the ion concentra- 
tion, which we measured to be 0.156 ± 0.008 kg/Mol for 
CuS04 using a density measurement instrument DMA 
5000 from Anton Paar. /j,a and represent the mobil- 
ity of the anions (8.3 10"^ m^/sV) and cations (5.6 10"® 
m^/sV), respectively. F is the Faraday constant (9.6 10''' 
As/Mol), g the acceleration due to gravity and Zc the 
charge number of the cation (2). represents the dy- 
namic viscosity (10~^ kg/ms) of the solution and D the 
ambipolar diffusion constant (8.6 10~^° m-^/s for CUSO4 
p6| ). Because D and are weakly concentration depen- 
dent, we assign them errors of 10 % in the subsequent 
calculations, further on we assume a 5 % error in the de- 
termination of d. Eq. H assumes, that the glass plates are 
at z = ztd/2 and y — at the anode. 

Within this theory the extension L of the convection 
roll is given by the point, where Vy{y, t) drops to zero: 



L = k2Vt 



(8) 



The maximal velocity i^max occurs at the heights z = 
±d/2-\/3, in the immediate neighborhood of the electrode 
(2/ = 0), and is time independent: 



= 0.63 



1/3 



Zc F (^a + Mc) V 



(9) 



However the authors state that due to problems with the 
boundary conditions at the electrode their solution might 
be not applicable in the region < y < d. 



2. Test of the theory 

In principle the theory should describe directly our ex- 
perimental results, however due to the uncertainties in 
the parameters we decided to perform a fit. Because we 
measure a two-dimensional projection of velocities, we 
concentrate on an average velocity Wavg(y,i)j which we 
calculate from all particles with absolute velocities of at 
least half the velocity of the fastest particle in this dis- 
tance. In Fig |l^ this corresponds to all particles not 
lying between the two solid lines. This restriction is nec- 
essary to remove the contribution of particles which have 
already settled to the bottom glass plate. 
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From Eq. |^ we derive our fit function: 



l) 12 V i 



(10) 



The fit of Eq. |o| to Wavg(y,0 is successful for the whole 
run of the experiment, especially the roll length L is de- 
termined very precisely. Fig. |ll| illustrates the temporal 
evolution of L obtained from the fit. After 70 s i^cexp be- 
comes constant at a value of 24.7 ± 0.5 fim/s. Integrating 
Eq. H yields wo.thco = 30.8 ±2.1 /im/s. 
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FIG. 11: Length L of the convection roll at the anode. L 
was determined by fitting «avg(j/, t) with Eq. The two fits 
to the left and right of the dashed line have slopes of 0.7 ± 
0.01 and 0.543 ± 0.001 where the errors are the asymptotic 
standard errors derived from the covariance matrix of the fit 
[ p?! . The experiment is the same as presented in Fig. [ic| . 



The next step is testing Eq. ^, which predicts the 
growth law of L. A fit of our experimental L for times 
greater 75 s using equation L{t) = k2jiti^ yields a slope 
of 0.543 ± 0.001, which is slightly above the square root 
law suggested by Eq. ^. Higher exponents have also been 
reported by other groups: Argoul et al. found 0.56 
± 0.01 and Dengra et al. [|| measured 0.54 ± 0.02. The 
coefficient fejit was found to be 141 ± 1 um/s^/^. If we 
insert our experimental parameters in Eq. |7| we obtain k2 
= 134 ± 12 ^m/s^/^ in agreement with the fit. 

Fig. ^ (b) shows the experimental results for the max- 
imal fluid velocity, which is located in the vicinity of the 
anode. After a sharp rise at the beginning of the experi- 
ment follows a slightly inclined plateau. A linear fit yields 
a velocity of 25.5 ±0.1 /im/s for i = 0, which increases 
at a rate of 1 % per minute. Inserting the experimental 
parameters into Eq. ^ results in a constant velocity of 
37.1 ± 2.5 /im/s. This discrepancy can be attributed to 
the unphysical boundary conditions used in the model. 

Summarizing it can be stated, that the theory 
presented in [ p^ provides a qualitatively and semi- 
quantitatively good description of the anodic convection 
roll in the DHS regime. 



3. Initial phase of development 

Apart from the DHS regime. Fig. |ll| also shows a 
growth law L ^ t^-'' for times between 12 s and 75 s. 
The exponent of 0.7 ± 0.01 indicates a faster growth 
originating from a different mechanism. In the very be- 
ginning of an experiment convective fluid transport is 
faster than the diffusive equalization between the copper 
ion enriched electrolyte and the bulk electrolyte. There- 
fore the concentrated electrolyte at the electrode sinks 
down and spreads along the bottom plate without sig- 
nificant mixing. In this so called immiscible fluid (IF) 
regime the length L of the convection zone is expected to 
grow with t'^/"" . This has been shown for gravity currents 
psi , E9| and was successfully adapted for electrodeposi- 
tion^ |l^. As explained in detail below, the range of 
applicability of this theory requires in our case L < 470 
/xm, which is out of our measured range. Thus we ob- 
serve a transitional period between the IF and the DHS 
regime and not the IF regime itself. 

To distinguish between the IF and the DHS regime 
a scaling analysis based on the vertical diffusion time 
was used in [^. We would like to advocate a differ- 
ent approach, using the similarity of the driving mech- 
anism with the well-investigated case of a side heated 
box filled with fluid [|o[ |3l|, ||. While in this case the 
density changes are due to thermal expansion, there ex- 
ist also two flow regimes: In the so called convective 
regime small layers of fluid spread along the confining 
plates with a stagnant core in the middle of the cell. The 
conductive regime is characterized by a cell filling convec- 
tion roll and iso-density lines which are almost vertical. 
These flows can be described by the Rayleigh number 
Ra = aAT pgS / KT] and the aspect ratio A = d/W . a 
is the thermal expansion and k the thermal diffusivity of 
the fluid. W denotes the distance, AT the temperature 
difference between hot and cold side wall. Boehrer |3^ 
pointed out, that RaA^, which equals the ratio between 
the timescales for vertical diffusion and horizontal con- 
vection, is the dimensionless control parameter of this 
transition. For high values of RaA^ one observes the 
convective regime, for small values the conductive one. 

If we transfer this analysis to our situation, we have 
to substitute the thermal density difference aATp with 
the density difference due to concentration changes 
and the thermal diffusivity k with the ambipolar diffu- 
sivity D. This yields a concentration-dependent Rayleigh 
number Ruc expressed by: 



Rac = 



dp 

dc 



Co .9 



d3 



D 1] 



(11) 



The aspect ratio is calculated using the length of the con- 
vection roll: A = d/L. Analyzing the results presented in 
Fig 10 of Ref. |9| this interpretation provides a necessary 
condition given by: 



Rac A^ > 1000 



(12) 
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to observe the IF regime. In our experiment Ruc will 
only be larger than 1000 for L < 470 /im , which is out 
of our measured range. 



B. Cathode Results 

At the cathode the situation is more complex due to 
the growing deposit. In Fig. ^ the soHd line describes 
the position of the most advanced point of the deposit 
in a width of 2.3 mm in x-direction, while the dotted 
line corresponds to the minimum in the same interval. 
The border of the convection roll was measured using a 
threshold: the open circles mark the foremost position, 
where the Vy component of a particle exceeded 4.3 //m/s. 

While the convection roll develops immediately, no 
growth of the deposit is observable in the first 40 s, be- 
cause the copper deposits in a planar compact way, which 
is not observable with our optical resolution. During this 
so called Sand's time, the ion concentration at the cath- 
ode drops to zero, which subsequently destabilizes the 
planar growth mode pO[ |. 

In the next phase (40 s < t < 280 s) a depth of the 
deposit (distance between the most advanced and most 
retarded parts of the growing deposit) becomes measur- 
able and finally reaches a constant size. The advancing 
deposit significantly compresses the size L of the convec- 
tion roll as given by the distance between the dots and 
the solid line. Within this time interval the hydrogel layer 
is established, which can be seen by visual inspection. 

The third phase {t > 280 s) is characterized by the 
appearance of the finger development. The front mini- 
mum (dashed line) and maximum (solid line) in Fig. |l^ 
coincide to the finger tip and the neighboring valley. The 
length L of the convection roll in front of the finger tip 
tends to converge to a constant size, which has also been 
observed in the absence of a gel ||9|, p2[ . The disper- 
sion relations studied in chapter piBf were obtained in 
this phase. 

Fig. |l^ is devoted to a comparison between the flow be- 
havior at the two electrodes. The developments of L dis- 
played in Fig. |l^ (a) differ significantly from each other. 
However the comparison of the maximal fluid velocities 
in Fig. |l^ (b) show similarities with respect to the abso- 
lute value and the approximate temporal Constance. The 
fit at the cathode yields a velocity of 23.6 ± 0.3 fim/s for 
t = s, which decreases 0.9 % per minute. 

From Fig. |l3| (a) we infer, that the theoretical descrip- 
tion presented in | [l9t cannot be applied to the cath- 
ode. Indeed three prerequisites of the theory are not ful- 
filled: Most importantly, it does not consider the moving 
boundary originating from the growth process. Moreover 
this model can lead to unphysical negative concentrations 
at the cathode due to its inherent simplifications. Finally 
the theory is two-dimensional in the y-z plane and there- 
fore not able to describe the influence of the ramified 
deposit, which evolves in the x-y plane. 

Remarkably enough our measurements show that the 
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FIG. 12: Temporal evolution at the cathode. The solid line 
describes the position of a finger tip, the dashed line of a 
neighboring valley. The circles mark the leading edge of the 
convection roll. The experimental conditions are: U — 12 Y 
and d = 300 /im. 
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FIG. 13: Comparison of the flow behavior at the two elec- 
trodes under identical experimental conditions. Filled circles 
correspond to the anode (Fig. ^l|), open circles to the cathode 
(Fig. y_2|). Part (a) shows the temporal evolution of the overall 
length L of the convection roll. Part (b) gives the temporal 
evolution of the maximal velocity in y direction inside the 
rolls. The solid lines are fits to all data points with t > 100 s. 
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presence of hydrogel does not prevent convection. For a 
better understanding of the contribution of the flow field 
to the finger morphogenesis, we visualize it in Fig. |lj for 
two different times of the same experiment. The solid 
line denotes the interface of the deposit, and the arrows 
indicate the velocity of individual particles. It is appar- 
ent that convection is restricted to a small zone in front of 
the growing deposit. The hydrogel occurs also in the im- 
mediate vicinity of the front, however its extension could 
not be investigated in any detail here. 
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FIG. 14: Flow field at (a) 493 s and (b) 760 s after the begin- 
ning of the experiment. All particles detected in one image 
with a velocity > 4.8 jim/s were included. Arrowheads have 
been omitted for clarity, the length of the lines corresponds 
to the way the particles travel in 6 s. {U = 12 V, d — 300 



roll and is displayed as the thick line in Fig. 

Fig. |l^ (b) reveals that the convection field has a mis- 
match of about 15 degrees to the front during the initial 
phase of finger development. Thus we conclude that the 
development of the convection field lags behind the devel- 
opment of the front. After the fingers are fully developed, 
the convection field readjusts again perpendicularly to 
the interface as shown in Fig. |l^ (c), which is a sign of 
the concentration gradient adapting to the geometry of 
the deposit. 

Due to the lack of comparable measurements of other 
electrodeposition systems we can not judge if this eflFect 
is due to hydrogel or a generic feature of the cathodic 
convection roll. 
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To examine the orientation of the flow field with re- 
spect to the interface we first computed the distribution 
of the angle a between the velocity vectors Vj of the par- 
ticles and the y axis. These data correspond to the thin 
lines shown in Fig. ^ The distribution clearly broadens 
with time and exhibits two distinct maxima for the fully 
developed fingers analyzed in Fig. |l^ (c). 

Then we identified for each particle an associated point 
of the finger envelope by prolonging v^. The angle be- 
tween the normal vector at this meeting point and the y 
axis is labeled (3. The angle difference a — /? is a measure 
for the mismatch between the deposit and the convection 



FIG. 15: Angle histograms for the times (a) 3-97 s, (b) 493- 
651 s and (c) 730-918 s. The thin line gives the distribution 
of the angle a between the velocity vector v of a particle and 
the y axis. Prolonging v we determine the meeting point with 
the finger surface. Its normal vector encloses the angle /3 with 
the y axis. The thick line shows the distribution of the angle 
difference a — j3. These data correspond to the experiment 
presented in Fig. O. 
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V. SUMMARY AND CONCLUSIONS 

We measured the dispersion relation in the hnear 
regime of the finger morphology and their dependence 
on cell thickness and applied potential. By means of 
a textured electrode, we were able to measure negative 
growth rates. The striking feature of the smooth fin- 
ger envelope is connected with the existence of a lim- 
ited band of wavenumber between an fccrit with posi- 
tive growth rates. The damping of all perturbations with 
higher wavenumbers can be attributed to an effective sur- 
face tension associated with a hydrogel boundary in front 
of the deposit A fit of the dispersion relation yields some 
estimates for the effective surface tension and the viscos- 
ity of the copper hydrogel. 

Furthermore we performed PIV measurement at both 
electrodes. At the anode we could confirm the i^^^ growth 
law for the length L of the convection roll. We deter- 
mined L by fitting the suggested analytic expression for 
the velocity field and could therefore successfully test the 
model proposed by Chazalviel and coworkers. 

At the cathode the maximal fluid velocity inside the 
convection roll is of the same order as at the anode, but 
the temporal evolution of L differs strongly. An analy- 
sis of the orientation of the velocity vectors reveals the 
existence of some mismatch between the development of 
the convection roll and the deposit front while the sys- 
tem is in the linear regime. In the fully developed finger 
regime, the velocity vectors are again perpendicular to 
the envelope of the deposit. 

Our results provide a reasonable explanation why in 
the absence of gravity-driven convection rolls the finger- 



ing instability cannot be observed: without convective 
mixing, the concentration gradient at the hydrogel inter- 
face can be assumed to be steeper, which will increases 
the effective surface tension. In consequence fccrit shifts 
to smaller wavenumbers and the overall growth rates de- 
crease, which will suppress the evolution of fingers. 

An alternative explanation assumes two zones within 
the hydrogel layer. One part of the hydrogel in immediate 
vicinity of the deposit will be mixed by the convection 
roll and the consequential shear thinning will decrease 
its viscosity. In front of it there is a zone of quiescent 
hydrogel of higher viscosity, at the interface between this 
two the instability takes place. In this scenario fccrit is 
determined by the length of the convection roll L. 

A fully quantitative theoretical analysis remains to be 
done. 
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